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Abstract 

We compare the two-loop corrections to the finite part of the light-like hexagon Wilson loop 
with the recent numerical results for the finite part of the MHV six-gluon amplitude in M = 4 
SYM theory by Bern, Dixon, Kosower, Roiban, Spradlin, Vergu and Volovich (]arXiv:0 803. 1465 
[hep-th]) and demonstrate that they coincide within the error bars and, at the same time, they 
differ from the BDS ansatz by a non-trivial function of (dual) conformal kinematical invariants. 
This provides strong evidence that the Wilson loop / scattering amplitude duality holds in planar 
M = 4 SYM theory to all loops for an arbitrary number of external particles. 
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1 Introduction 



In an important recent development in the study of the AdS/CFT correspondence, Alday and 
Maldacena proposed [TJ the strong coupling description of planar gluon scattering amplitudes 
in the maximally supersymmetric TV = 4 Yang-Mills (SYM) theory and were able to make a 
direct comparison with a prediction based on weak coupling results for the same amplitudes. At 
weak coupling, following years of intensive studies of gluon scattering amplitudes [21 [31 0] , the 
conjecture was put forward by Bern, Dixon and Smirnov [5] (see also an earlier version in [6]) 
that the maximally helicity-violating (MHV) planar gluon amplitudes in M = 4 SYM have a 
remarkably simple all-loop iterative structure. In general, these amplitudes have the following 
form: 

lnA4 MHV ) = [IR divergences] + F^^fa, ...,p n ;a) + 0(e) . (1) 

Here .Mi MHV ^ is the color-ordered planar gluon amplitude, divided by the tree amplitude. The 
first term on the right-hand side describes the infrared (IR) divergences and the second term 
is the finite contribution dependent on the gluon momenta pi and on the 't Hooft coupling 
a = g 2 N/(87i 2 ). The structure of IR divergences is well understood in any gauge theory [TJ. In 
particular, in theories with a vanishing beta function like M = 4 SYM, the leading IR singularity 
in dimensional regularization is a double pole, whose coefficient is the universal cusp anomalous 
dimension appearing in many physical processes [HI El QUI ED] • The BDS conjecture provides an 
explicit expression for the finite part, i^ MHV ^ = F^ DS \ for an arbitrary number n of external 
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gluons, to all orders in the coupling a. Remarkably, the dependence of Fn on the kinematical 
invariants is described by a function which is coupling independent and, therefore, can be deter- 
mined at one loop. At present, the BDS conjecture has been tested up to three loops for n = 4 
[2] and up to two loops for n = 5 [T2]. The explicit investigation of the conjecture for n = 6 at 
two loops is the subject of this paper and of the parallel paper [T5] . 

As mentioned earlier, after the work of Alday and Maldacena it became possible to test the 
BDS conjecture at strong coupling. According to their proposal, at strong coupling the planar 
gluon amplitude is related to the area of a minimal surface in AdSs space attached to a specific 
closed contour C n , made out of n light-like segments defined by the gluon momenta 

x i ~ x i+i = Pi (with the cyclicly condition x n+1 = x\), 



Jq 2 N 

\nM n = -^l—A min (C n ). (2) 

For n = 4 the minimal surface A min (C4) was found explicitly in [Tj, by making use of the conformal 
symmetry of the problem. With the appropriate AdS equivalent of dimensional regularization, 
the divergent part of lnAdi has the expected pole structure, with the coefficient in front of the 
double pole given by the known strong coupling value of the cusp anomalous dimension. Most 
importantly, the finite part of In A^4 is in perfect agreement with i^ BDS ) f r0 m the BDS ansatz. 
For n > 5 the practical evaluation of the solution of the classical string equations turns out to be 
difficult, but it simplifies significantly for n large [TJ]. In the limit n — > oo the strong coupling 
prediction for \nAi n disagrees with the BDS ansatz. This indicates [TJ] that the BDS conjecture 
should fail for a sufficiently large number of gluons and/or at sufficiently high loop level. 

Alday and Maldacena pointed out [TJ that their prescription (T5]) is mathematically equivalent 
to the strong coupling calculation of the expectation value of a Wilson loop W(C n ), defined on 
the light-like contour C n [TBI, [TB] . This should not come as a total surprise, since the intimate 
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relationship between the infrared divergences of the scattering of massless particles and the 
ultraviolet divergences of Wilson loops with cusps is well known in QCD [HI HD]. Inspired 
by this, in [17] three of us conjectured that a similar duality relation between planar gluon 
amplitudes and light-like Wilson loops also exists at weak coupling. We illustrated this duality 
by an explicit one-loop calculation in the simplest case n = 4. This was later extended to the 
case of arbitrary n at one loop in [18j. If we write the log of the Wilson loop in the following 
way, 

In W{C) = [UV divergences] + F^ t) {x 1 , ...,x n ;a) + 0(e) , (3) 

then the duality relation identifies the finite parts of the two objects up to an additive constant, 
once one imposes the relations Vi =x i~ x i+n 

F( MHV )=Fi WL )+ const. (4) 

This property is extremely non-trivial. 

Further evidence in favor of the duality relation (jlj) at weak coupling came from two-loop 
calculations of lnW(C n ) for n = 4 [20] and n = 5 [21 J . Our results were in agreement with 
the two- loop MHV gluon amplitude calculations [6j [12], and hence with the BDS ansatz for 
n = 4,5. Furthermore, in [2U] we proposed and in [2T] we proved a conformal Ward identity for 
the light-like Wilson loop W(C n ), valid to all orders in the coupling. It fixes the functional form 
of the finite part of hxW(C n ) for n = 4 and n = 5, up to an additive constant, to agree with the 
conjectured BDS form for the corresponding gluon amplitudes^ However, for n > 6, although the 
BDS ansatz Fn BDS ^ satisfies the conformal Ward identity, F„ L ^ is allowed to differ from _F 1 i BDS ' ) 
by an arbitrary function of conformal invariants (for n = 6 there are three such invariants). This 
result provided a possible explanation of the BDS conjecture for n = 4, 5 (assuming that the 
MHV amplitudes have the same conformal properties as the Wilson loop, see the discussion in 
Sect. 2.4 below), but left the door open for potential deviations from it for n > 6. To verify 
whether the BDS conjecture and/or the proposed duality relation (J4j) still hold for n = 6 to 
two loops, it was necessary to perform explicit two-loop calculations of the finite parts of the 
six-gluon amplitude i^ MHV ^ and of the hexagon Wilson loop Fq WL \ 

In a recent paper [24] we reported on the two-loop calculation of Fg WL ' 1 and found that it 
differs from the BDS ansatz, 

F (WL) = F (BDS) + Rq (5) 

Here, in complete agreement with the conformal Ward identity, Rq is a non-trivial 'remainder' 
function of three conformally invariant combinations of the kinematical variables. As was empha- 
sized in [24J, were the duality relation (T4J) to hold for n = 6, the function R 6 would describe the 
discrepancy between the BDS ansatz and the scattering amplitude. A parallel two-loop six-gluon 
amplitude calculation has been undertaken by Bern, Dixon, Kosower, Roiban, Spradlin, Vergu 
and Volovich and the results became available very recently [13J. 

The detailed numerical comparisons of the two calculations, described in this paper and in 
the parallel publication |13j, shows that, firstly, the BDS ansatz fails for n = 6 at two loops and, 
secondly, the duality with Wilson loops is preserved, 

F (MHV) = F (WL) _ Ce(Q) > (6) 



1 Later, similar Ward identities were also obtained at strong coupling using the AdS/CFT correspondence in 

Refs. 122 [231 . 
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where cg(a) is a constant. We consider this as very strong evidence that the duality relation (jlj) 
should hold for arbitrary n to all orders in the coupling. 

The real challenge now is, in our opinion, to find out the deep reason behind the surprising 
duality (jlj) between two apparently unrelated objects in the Af = 4 theory. One might speculate 
that the scattering amplitudes and the light-like Wilson loops share the same (probably infinite) 
set of symmetries, of which (dual) conformal symmetry is just the most visible part. If so, these 
symmetries may completely fix the form of the quantities on both sides of the duality relation. 
This could be the manifestation of some new type of integrability of Af = 4 SYM theory. We 
would like to stress that the duality relation (jlj) only holds in the planar limit. Indeed, the 
known two-loop non-planar contributions to the four-gluon amplitude [3] appear to break dual 
conformal symmetry, which, as we have shown, is an important ingredient of the duality. Also, 
recent studies in the AdS/CFT correspondence suggest [221 [2S1 [2Z| that at strong coupling the 
duality relation (J3j) can be extended to scattering amplitudes involving gluinos (both in the 
M = 4 SYM and its deformed versions), as well as to the on-shell matrix elements of conserved 
currents [14J. It would be interesting to verify whether the same relation holds at weak coupling. 
This question is beyond the scope of the present paper and deserves further investigation. 

We would like to point out that a weaker form of the duality (jlj) has already been observed 
in QCD in the special, high-energy (Regge) limit s ^> — t > for the four-gluon amplitude up to 
two loops [10]. The same relationship holds in any gauge theory ranging from QCD to Af = 4 
SYM. The essential difference between these theories is that in the former case the duality is only 
valid in the Regge limit, whereas in the latter case it is exact in general kinematics. Moreover, 
in Ref. [T7], based on the BDS ansatz, three of us argued that the four-gluon amplitude is Regge 
exact in Af = 4 SYM. Namely, the contribution of the gluon Regge trajectory to In A4^ coincides 
with its exact expression evaluated for arbitrary values of s and t, up to terms vanishing as 
e — > 0. This property allowed us to obtain the explicit expression for the three-loop gluon Regge 
trajectory in Af = 4 SYM [17J. The two-loop correction to this trajectory was found to be 
in agreement with the results of Ref. |28j and the three-loop correction was later confirmed in 
Refs. [23 1301 ED E2]- We would like to mention that a thorough analysis of the Regge limit of 
planar multi-gluon amplitudes in Af = 4 SYM was recently performed in Refs. [301 EI] arid it 
provided further evidence that the BDS ansatz needs to be corrected [3T] . 

The paper is organized as follows. In Section [2] we explain in detail the proposed duality 
between planar gluon amplitudes and light-like Wilson loops. We summarize some basic facts 
about planar gluon amplitudes, focusing on the structure of their infrared divergences and on 
the BDS ansatz for the finite part. Next we describe the light-like Wilson loop, the structure of 
its ultraviolet singularities and the anomalous conformal Ward identity for its finite part. Then 
we state the duality relation and discuss its consequences for the six-gluon amplitude. Section 
3 contains a description of our two-loop calculation of the hexagon Wilson loop W(Cq) (already 
reported in [24]). We separate the divergent and finite part of individual Feynman diagrams by 
employing the 'subtraction procedure' proposed in [21], and then demonstrate that the divergent 
part of In W(C 6 ) is of the expected form. We work out the representation for the finite part i^ WL ^ 
in the form of convergent multiple parameter integrals which can easily be evaluated numerically 
for given kinematical configurations. In Section H] we perform a detailed numerical comparison 
of our results of the hexagon Wilson loop calculation with the BDS ansatz for n = 6 and with 
the results of the parallel six-gluon calculation [TJ]. Section 5 contains concluding remarks. 
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Figure 1: The conjectured duality relation between the gluon scattering amplitude AA n and the Wilson 
loop W(C n ). The dashed lines depict gluons and the double lines the integration contour Cq. The 
momenta of the incoming gluons are identified as the light-like segments of the integration contour, 
Pi ~ Xi - x i+ \. 

2 Planar gluon amplitude/ Wilson loop duality 

Recent studies revealed that the gluon scattering amplitudes have a number of remarkable prop- 
erties in the M = 4 SYM theory. To describe them, we first recall some general features of such 
amplitudes. 



2.1 Planar amplitudes 

In a generic Yang-Mills theory with an SU(N) gauge group, the scattering amplitude of n gluons 
can be decomposed into color-ordered partial amplitudes multiplied by the corresponding color 
structure [33]. In the planar limit, the dominant contribution only comes from the single-trace 
color structures leading to 

A n ({ Pu h u a t }) = 2 n ^g n ~ 2 tr[t a ^K.A a ^]A n (a(l h \...,n hn ))+0(l/N 2 ), (7) 

where each gluon is characterized by its on-shell momentum (p? = 0), helicity hi = ±1 
and color index Oj. Here the sum runs over all possible non-cyclic permutations a of the set 
{1, . . . ,n} and the color trace involves the generators t a of SU(N) in the fundamental represen- 
tation normalized as tr(t a t b ) = \5 ab . All gluons are treated as incoming, so that the momentum 
conservation takes the form Y17=i Pi = ®- 

The planar scattering amplitude ([7j) is uniquely determined by the set of color-ordered am- 
plitudes A n . It follows from the supersymmetric Ward identities [34] that the amplitude A n 
vanishes to all orders when either all the external gluons or all gluons but one have the same 
helicity. When two gluons have, e.g., negative helicity and the remaining n — 2 gluons have 
positive helicity, the so-called maximally helicity- violating (MHV) amplitudes have the following 
remarkably simple form in the spinor helicity formalism: 

A (1 + m- m~ n + ) - • (m 1 m 2 ) A u(mhv) /o\ 
An[L ••• mi ••• m2 ) ~\l2)(23)...(nl) Mn • W 
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Here the helicity structure is described by massless Weyl spinors u±(p) of momentum p and chiral- 
ity ±1; they form Lorentz invariant 'inner products' of the type (j k) = (j~\k + ) = u^(pj)u + (pk). 
The Lorentz scalar Jvt n = 1 + aM. n -\ + 0(a 2 ) does not depend on the positions mi j2 of the 
negative- helicity gluons [35J. According to (|7j) and ©, the perturbative corrections to the planar 
MHV amplitude are determined to all loops by the single function A^n MHV ^ of the Mandelstam 
variables. This function (more precisely, its logarithm) is the main object of interest in this 
paper. 

It is worthwhile mentioning that for scattering amplitudes with more than two negative 
helicities (next-to-MHV amplitudes and their generalizations) the situation is more complex. 
At tree level, in the spinor helicity formalism these amplitudes are expressed in terms of various 
Lorentz structures built from spinors [33]. At one-loop level, however, the perturbative corrections 
to the non-MHV amplitudes induce new Lorentz structures [1]. As a result, in distinction with 
(jSD, the non-MHV amplitudes do not admit a simple factorized form. This is the reason why the 
duality with Wilson loops that we discuss here only concerns MHV amplitudes. 



2.2 Planar MHV amplitudes 

The gluon scattering amplitudes suffer from infrared divergences. In supersymmetric gauge 
theories, they can be regularized using the dimensional reduction scheme (DRED) with D = 
4 — 2eiR and em < 0. @ In this case, the IR divergences appear as poles in eiR. At a given loop 
order /, the maximal order of the poles is 21. In a generic Yang-Mills theory, the IR divergent 
part of the planar scattering amplitudes has a universal form which is sensitive neither to the 
helicities of the scattered particles, nor to their type (gluons, fermions, scalars). More precisely, 
the IR divergences exponentiate in the all-loop planar amplitude and, as a consequence, they can 
be factorized into a universal divergent factor. In application to the MHV planar amplitudes in 
M = 4 SYM theory, this property allows one to decompose In .Mi MHV ^ into divergent and finite 
parts as follows, 

ln _ M (MHV) =Zn + F (MHV) +0(effi); (g) 

In a theory with a vanishing beta function, the IR divergent part Z n takes the particularly simple 
form 



4£' 



cusp 



Q(l) \ »* / t [2] ^ ~ klR 
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where tj 2 ' = s^+i = (pi + Pi+i) 2 is the invariant mass of two adjacent gluons with indices % and 
? + l subject to the periodicity condition i + n = i. The IR cut-off /z 2 R is related to the dimensional 
regularization scale /i 2 as 

/4 = 47re^V, (11) 

where 7^ is the Euler constant. 

We did not put the superscript MHV on Z n in order to indicate that it has a universal form 
for all planar amplitudes in the M = 4 SYM theory. The coefficients of the IR poles in fflUj) are 
the expansion coefficients of the cusp anomalous dimension r cusp (a) = ^rlusp and of the 



2 More precisely, a modification of the DRED scheme, the so-called four-dimensional helicity (FDH) scheme, 
has been used in [36l. 
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so-called collinear anomalous dimension G(a) 
in the DRED scheme by 



Ylili a l G( l \ To two- loop order, they are given 



cusp (o) 



2a - 2( 2 a 2 + 0(a 3 



G{a) = -C 3 a 2 + 0(a 3 



(12) 



By definition [37J El [38], the cusp anomalous dimension T 



cusp V 



describes specific ultraviolet 
divergences of a Wilson loop evaluated over a contour with a cusp. Its appearance in the infrared 
divergent part of the scattering amplitude (fTO]) is not accidental, of course. It has its roots in the 
deep relation between scattering amplitudes in gauge theory and Wilson loops evaluated over 
specific contours in Minkowski space-time, defined by the particle momenta [HEME]]. It should 
be mentioned that this relation is not specific to M = 4 SYM and it holds in any gauge theory, 
including QCD. The two- loop expression for r cusp (a) in a generic (supersymmetric) Yang- Mills 
theory was found in Refs. [351 11]. We remark that in M = 4 SYM theory, r cusp (a) is known 
at weak coupling to four loops [5H], and there is a conjecture for it to all loops [40J. At strong 
coupling, the solution to the BES equation proposed in [ID] produces a strong coupling expansion 



of r cusp (a) [JTJ H21 H3j. The first few terms of this expansion are in agreement with the existing 



quantum superstring calculation of Refs. [HJ H3 H6]. The non- universal collinear anomalous 
dimension G(a) is known to four loops at weak coupling SZ|. 

Let us now examine the finite part of the MHV planar amplitude Fi MHV ^ . In a generic Yang- 
Mills theory, it has a much more complicated form compared to the divergent part. Surprisingly, 
this finite part becomes considerably simpler in M = 4 SYM. Previously, the finite parts of the 
following MHV planar amplitudes have been calculated: 

• n — 4 gluon amplitude up to three loops [5J, [5] ; 

• n — 5 gluon amplitude up to two loops [12J; 

• n > 6 gluon amplitude at one loop [2]. 

These calculations revealed a remarkable iterative structure of F,( MHV \ They led to the formu- 
lation of the Bern-Dixon-Smirnov ansatz (BDS) which provides a conjectured expression for the 
finite part of the MHV planar amplitudes valid for an arbitrary number of gluons n > 4 and to 
all orders in the coupling. The BDS ansatz reads 



p(BDS) _p 

2 



cusp (a) T n + c(a) 



(13) 



where r cusp (a) 



is the cusp anomalous dimension defined in f[T2l and c(a) = — \(,\a 2 + 0(a 3 ) is 
a universal (independent of n) additive constant. The dependence on the kinematical invariants 
is described by the function T n . According to the BDS conjecture, this function is coupling 
independent, and can thus be determined at one loop. For our purposes here we only need its 
explicit expressions for n = 4, 5, 6: 



^ 4 =-l n -(-l r )+4C 2) 



(14) 



Ft 



1 5 

i=l 

i 6 r 



J2] 



[2] 



i+3 



In 



J2] 



v i+2 



In 



J2] 

'i+1 

+ [3] 



3/ 
+ 2 C2 



4 



t\ 



.13] 
b i+l 



Lio 

2 



J2U2] 
J3U3] 



4 C2 



(15) 



(16) 
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where t[ = (pi + ... + pi +r -i) 2 are the Mandelstam kinematical invariants. 
2.3 Light-like Wilson loops 

Let us now turn to the description of the light-like Wilson loops, which are the counterparts of 
the MHV planar amplitudes in the duality relation discussed in this paper. In the M = 4 SYM 
theory with an SU (N) gauge group they are defined as 

W(Cn) = ^(0|TrPexp (i £ dx^xfj |0) , (17) 

where the gauge field A^x) is integrated along the contour C n = UILi ^ ma de out of n light-like 
segments joining the cusp points xf (with i = 1, 2, . . . , n) 

ii = {^(Ti) = nxf + (l - Ti)x>t +l \ r, t e [0, l]} , (18) 

such that the tangent vectors <9 Ti x M (r,;) = x^ i+1 are light-like, xf i+l = 0. The symbol P indicates 
the ordering of the SU(N) indices along the integration contour C n . 

The Wilson loop (TTTT) is a gauge invariant quantity depending on the integration contour C n . 
As a function of the cusp points, it is invariant under their cyclic permutations and flips, 

W(xt, x 2 , . . ■ , x n ) = W(x n , xi, . . . , x n _i) = W(x n , x n -i, . . . , xi) . (19) 

If the contour C n did not have cusps, the Wilson loop would be a finite quantity in D = 
4 [371 HBJ HE]. Characteristic feature of the light-like contour C n is that it maps to a similar 
contour under 5*0(2,4) conformal transformations. As a result, a finite Wilson loop would be 
a conformal invariant function of Xi [TJ [20]. We would like to stress that, contrary to the gluon 
scattering amplitudes, the Wilson loop fflTj) is defined in configuration space and the conformal 
transformations act on the four- vectors xf defining the positions of the cusp points in Minkowski 
space-time. Due to the presence of cusps on the integration contour C n , the Wilson loop ( 1T7I) 
has specific ultraviolet divergences [37J HHJ [HJ [381 EHJ M which make the conformal symmetry of 
W{C n ) anomalous. 

To regularize the cusp singularities, we use dimensional reduction with D = 4 — 2e uv and 
e uv > (notice the sign difference compared to e IR ). Like the scattering amplitude, the Wilson 
loop can be split into a divergent and a finite part, 

In W(C n ) = + Fi w ^ + 0(e uv ) . (20) 

The divergent part Z^ WL ' ) has the special form [19] 

4 i>i V €vv > L€vv J «=i 

where Tcusp are the expansion coefficients of the cusp anomalous dimension (Tl2|) and is the 
counterpart of the collinear anomalous dimension entering in (flUj) . 

T(a) =^^r (/) = -7C 3 a 2 + 0(a 3 ). (22) 
i>i 
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The UV cut-off /i uv is related to the dimensional regularization scale [i as (cf. (Hip) 

/4 v = vre^V- (23) 



The finite part of the Wilson loop, Fi WL \ does not depend on the renormalization scale \i 



UV 



and it is a dimensionless function of the distances with i, j = 1, . . . , n. Since the edges of C, 



are light-like, xf i+1 = 0, the non- vanishing distances are xfj with \i — j\ > 2. Despite the fact 
that the conformal invariance of W(C n ) is broken by the cusp singularities, it imposes severe 
constraints on Fi WL ' ) . We showed in [2TJI |2~T] that F^ L ^ has to satisfy the following anomalous 
conformal Ward identity: 

n i n / 2 \ 

^ F (WL) ^ (2*^ • ft - x^) Fi WL) = lr cusp (a) ^ In ( . (24) 

i=i / <=i V^-i.i+i/ 

The differential operator on the left-hand side is the conformal boost (special conformal trans- 
formation) generator . The right-hand side expresses the conformal anomaly due to the cusp 
singularities. As explained in [21] , it has a universal functional form, with the coupling depen- 
dence coming only through the cusp anomalous dimension r cusp (a). For n = 4 and n = 5 this 
relation is powerful enough to determine the all-loop expressions for F^ WL ^ and F^ WL ^ up to an 
additive coupling-dependent constant. The reason for this is that any potential solution of the 
homogeneous differential equation K^F^ = would be a conformal invariant. It is well known 
that such invariants take the form of cross-ratios x\^x\ x j [x\ k x 2 ^ . It is then immediately clear 
that one cannot build invariants from four or five points with light-like separations xf i+l = 0. 
However, they can be constructed starting from six points. In particular, for the hexagon Wilson 
loop W(Cq) there are three such cross-ratios, 

_ " L 13 X 46 „, _ X 24 x 15 „, _ X 35 X 26 /r,r\ 
U\ — — g 5-, «2 — — 2 o— , «3 — — q 2~ " V Z0 / 

x 14 x 36 x 25 x 14 x 36 x 25 

As a result, for n > 6 the general solution of the conformal Ward identity (124")) will contain an 
arbitrary function of the conformal cross-ratios. This function was calculated at two loops for 
n = 6 in [21]. In the present paper we give more details of this calculation and compare the 
result with the corresponding six-gluon MHV amplitude. 



2.4 Duality relation 

In this subsection we formulate and discuss the main point of the present paper - the proposed 
duality relation between the MHV planar amplitudes .A/fi MHV ^ and the light-like Wilson loops 
W(C n ). 

The conjectured duality states that in the planar M = 4 SYM theory the finite parts of the 
logarithms of the gluon amplitude and of the Wilson loop are equal (up to an inessential additive 
constant): 

f (MHV) =F (WL) + congt) (26) 

upon the formal identification of the external on-shell gluon momenta in the amplitude with the 
light-like segments forming the closed polygon C n (the contour of the Wilson loop), 

(27) 



8 



Thus, the Mandelstam variables for the scattering amplitudes tj = + . ■ .+pi + j~i) 2 are related 
to the distances xfj between two cusp points on the integration contour of W(C n ) as follows, 

W^fc := X 1,i+jl x \,k+i ■ (28) 

The divergent parts of the scattering amplitudes and the light-like Wilson loops are also 
related to each other but the relationship is more subtle since the two objects are defined in two 
different schemes (infrared regularization for the amplitudes and ultraviolet regularization for the 
Wilson loops), both based on dimensional regularization. From the discussion in Subsections 12.21 
and 12.31 we know that the leading IR divergence of the amplitude (the coefficient of the double 
pole in Eq. (TTUI) ) coincides with the leading UV divergence of the Wilson loop (the coefficient 
of the double pole e^y in Eq. ff2TT) ). since both are controlled by the universal cusp anomalous 
dimension r cusp (a). One can also achieve the matching of the coefficients of the sub leading simple 
poles corresponding to the (non- universal) collinear anomalous dimensions G(a) and T(a) from 
Eqs. dTJ) and J22j), respectively!! To this end one relates the parameters of the two different 
renormalization schemes as follows, 

4 +2 /4v := 4*7/4 e 7(a) , e uv := -e m e^ . (29) 

Here the functions 7(a) and 5(a) are chosen in a way to compensate the mismatch between G(a) 
and T(a), without creating extra /i-dependent finite terms. It is easy to check that these functions 
are the solutions to the equations 

7(«)r cU s P (a) + 5(a)f cusp (a) + G(a) + T(a) = , 

7(a)f cusp (a) + 2<5(a)f cusp (a) + G(a) +f(a) = 0, (30) 

where f cusp (a) = f Q a ^-T cusp (a') and f cusp (a) = J Q a cusp (a'), and similarly for G(a) and f (a). It 
should be stressed that this procedure is not analogous to comparing two different renormalization 
schemes for the computation of the same divergent object. It is rather a change of variables 
(regularization parameters) which allows us to compare two different objects computed in two 
different schemes. Another such change of variables is the identification (1271) of the particle 
momenta with the light-like segments on the Wilson loop contour. 

The duality relation (1261) was inspired by the prescription of Alday and Maldacena for comput- 
ing gluon scattering amplitudes at strong coupling [TJ] , which essentially recasts the amplitudes 
into light-like Wilson loops in the dual variables ff27|) . A priori, one would expect the strong 
coupling relation between gluon amplitudes and Wilson loops to receive 1 / a/A corrections, which 
might spoil the relation at weak coupling [H] . Nevertheless, in [17] three of us found that at one 
loop and for four points the Wilson loop and the gluon amplitude agree, which lead to the idea 
that the duality might also be true perturbatively. In [TS] the duality at one loop was shown 
to apply also to n-point amplitudes. Further evidence in favor of the duality was accumulated 
by the present authors by two-loop calculations of the four-point [20] and five-point [21] Wilson 
loops, which were shown to match the corresponding gluon amplitudes according to Eqs. (1261) 
and (|27j). 

3 We thank Paul Heslop for turning our attention to the incomplete discussion of this point in the first version 
of the present paper. We are also grateful to Lance Dixon for a discussion of the different physical interpretations 
of the IR and UV simple poles [39] • 
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What can we say about the proposed duality to all orders in the coupling a! Supposing the 
duality holds, one immediate consequence is that Fi MHV \ after making the change of variables 
(I27p . should satisfy the conformal Ward identity ([H]) . We stress the fact that, while the conformal 
properties of the Wilson loop are manifest, there is no obvious reason to expect that the matching 
gluon amplitudes should have any reasonable behavior under conformal transformations acting on 
the particle momenta. Yet, the BDS ansatz for the MHV amplitudes does have this unexpected 
conformal property. For instance, when rewritten in terms of the dual variables £j ( 1271) . the BDS 
ansatz for four (fT4"|) and five (fT5"i) gluons reads as follows: 

Ff DS) = 7 Wa) In 2 $f) + const , (31) 

Ff DS) = -k usp (a) X>(4^) In(^N^) + const . (32) 

It is easy to check that these formulae are indeed solutions of the Ward identity (I24|) . In fact, 
they are the unique solutions (up to an additive constant) for n = 4,5, so they give the all-order 
form for i^ WL ^ and F^ WL \ T n particular, this confirms the duality up to three loops for n = 4, 
since on the one hand, the BDS ansatz is known to be correct up to three loops in this case, and 
on the other hand, the n = 4 Wilson loop does satisfy the Ward identity to all orders. 

A hint at a possible source of this surprising conformal symmetry of the MHV amplitudes 
comes from the observation that all momentum integrals entering in the expression for .M4 up 
to four [5j [39] (and possibly even five [50]) loops are of the 'pseudo-conformal' (also referred to 
as 'dual conformal') type [5TJ [17] The conformal properties of such momentum integrals are 
revealed by rewriting them in terms of 'dual coordinates' according to ( J27I) . If their external legs 
are taken off shell, one can remove the dimensional regulator and the integrals become manifestly 
covariant under the action of the conformal group 5*0(2,4) on the dual coordinates (i.e., on the 
particle momenta). It should be made very clear that this symmetry of the integrals (broken on 
shell by the IR divergencies) is not related to the original conformal symmetry of the M = 4 SYM 
theory, acting on the gluon fields in the configuration space. It is also important to realize that 
this unexpected property of the momentum integrals contributing to the MHV gluon amplitudes 
does not automatically imply that the amplitude should satisfy the anomalous (dual) conformal 
Ward identity f[24"j) . This is only true for some very special combinations of such pseudo-conformal 
integrals (for example, all the integrals in the four-gluon amplitude up to four (five) loops appear 
with coefficients ±1 [391 EB EE3 E2]). The exact link between the pseudo-conformal property of 
the integrals and the dual conformal behavior of the amplitude is still unclear. 

At this stage one might suspect that the observed duality relation for n = 4, 5 is true only 
because both objects have the same conformal symmetry. To put it differently, the proposed 
duality might be reduced to the weaker (but still highly non-trivial) statement that the MHV 
gluon amplitudes have dual conformal symmetry. The first real test of the stronger form of the 
duality ( |26i) is provided by the case n = 6. 

4 The conformal properties of ladder (scalar box) multiloop integrals were first described by Broadhurst (see 
the first reference in [51]). 
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The BDS ansatz for the six-gluon amplitude ( 1T61) . rewritten in the dual variables Xj, reads 



ry 2 , y . 2 

1 In' 



i=1 L x "%i+3 x ^i,i+3 

1 . 

+ - In 



1/ T* ^ \ "1 / ,yj 2 ^y«2 

■ 2/ X «,«+3 \ J- T . x i,i+2 J/ i+3,i+5 



I V 'T*^ / O \ /-yj 2 ™ 2 

1 vx i+l,i+4 / ^ x x i,i+3 x i+2,i+5 



+ const . (33) 



Again, we can immediately verify that this is a solution of the Ward identity (1241) . However, 
unlike the cases n = 4, 5 where the solutions (I3TI) . (|32|) are unique up to an additive constant, 
this is no longer true for n > 6. As explained at the end of Subsection 12.31 for n > 6 the 
general solution of (1241) can contain an arbitrary function of the conformal cross-ratios (and of 
the coupling constant). Thus, for n = 6, we can expect a possible deviation from the BDS ansatz 
of the form 

F (WL) = ^(BDS) + R ^ ^ a) _ (34) 



The one-loop Wilson loop calculation of p8j has shown that the 'remainder' function Rq is 
just a constant at one loop, which is a confirmation of the Wilson loop/scattering amplitude 
duality going beyond the scope of conformal symmetry. However, one might suspect just a low 
loop-order 'accident'. The point is that the function R e must satisfy a further, rather powerful 
constraint, the so-called collinear limit [21 E] (see also Subsection 14.11 for further comments). 
It could be that due to the limited choice of loop integrals at this low perturbative level, the 
function (13"31 made of them is the only one satisfying both the conformal Ward identity (|24|) and 
the collinear limit. If so, at some higher perturbative level new functions with these properties 
might appear which could spoil the BDS ansatz and/or the Wilson loop/scattering amplitude 
duality. 

Indeed, our recent explicit n = 6 Wilson loop calculation [21] has shown that at two loops 
there exists a non-trivial 'remainder' function i? 6 satisfying both conditions. The crucial test 
then was to compare the results of our Wilson loop calculation with a parallel two-loop six-gluon 
amplitude calculation, in order to check whether the proposed duality between Wilson loops 
and gluon amplitudes continues to hold at this level (which, if true, automatically implies the 
breakdown of the BDS ansatz). The results of the six-gluon calculation have become available 
very recently [13]. As we demonstrate in Section HJ the detailed numerical comparisons of the 
two calculations shows that indeed the BDS ansatz fails at two loops whereas the duality with 
Wilson loops is preserved, 

F (MHV) = F (WL) _ C6(fl) (35) 

We consider this very strong evidence that the duality should hold to all orders in the coupling, 
although the 'remainder' function R n is likely to receive corrections at each loop order. 



3 Light-like hexagon Wilson loop 

The two-loop calculation of the light-like hexagon Wilson loop W(Cq) goes along the same lines 
as the analysis of the rectangular (n = 4) and pentagonal (n = 5) Wilson loops W(C n ) performed 
in Refs. [201 02], where the interested reader can find the details of the technique employed. 

As was explained in detail in [20J, the two-loop calculation of W(Cq) can be significantly 
simplified by making use of the non-Abelian exponentiation property of Wilson loops [53]. In 
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application to W(Cq), it can be formulated as follows: 

\nW(C 6 ) = £~ 2 C F w^ + (j^ C F Nw&> + 0(g") , (36) 

where C F = (N 2 — 1)/(2N) is the quadratic Casimir of the SU(N) in the fundamental repre- 
sentation and tin 1,2 ) are dimensionless (iV— independent) functions of the distances a; 2 - and UV 
cut-off fi 2 v . In calculating we do not rely on the planar limit. As a result, equation ( 1361) 
is exact in N. E We will take the planar limit later on, when we need to compare our result to 
the planar MHV amplitude. The one-loop correction was found in [18J. Its divergent part 
coincides with the one-loop contribution to Z n , Eq. (TlOl) . and its finite part coincides up to an 
additive constant with the BDS ansatz ffl3|) and satisfies the duality relation (j4j) at one loop. 

For finite N the two-loop corrections to W(Cq) involve two different color factors Cp and 
C F N. According to ( 1361) . the coefficient in front of g 4 C F /(47r 2 ) 2 in the two-loop expression for 
W(Cq) is given by (w^ 1 ') /2. Therefore, in order to determine the function one has to 
calculate the contribution to W(Cq) only from two-loop diagrams containing 'maximally non- 
Abelian' color factor C F N . This property significantly reduces the number of relevant two-loop 
diagrams. The function is gauge invariant, so in order to simplify the calculation we shall 
employ the Feynman gauge. In this gauge, some of the 'maximally non-Abelian' diagrams like 
those where both ends of a gluon are attached to the same light-like segment vanish by virtue of 
x "j j+i = 0- The corresponding Feynman diagrams have the same topology as for the rectangular 
(n = 4) and pentagon (n = 5) Wilson loops W(C n ) and they can easily be identified by applying 
the selection rules formulated in Ref. [20J. 

To summarize, in Fig. [2] we list all non-vanishing two-loop diagrams of different topologies 
contributing to w^ 2 \ The diagrams shown in Figs. [2](m) - (r) involve the three-gluon interaction 
vertex of the Af = 4 SYM Lagrangian. Their color factors equal C F N, and therefore they 
contribute to the function in (|36|) . The non-planar diagrams in Figs. GD^a) - (1) involve two 
free gluon propagators and their color factors equal C F (C F — N/2). To identify their contribution 
to w^ 2 \ we have to retain the maximally non-Abelian part only, that is, to replace their color 
factors by C F (C F —N/2) — > —C F N/2. Finally, the diagrams in Figs.[2](s) -(u) involve the one-loop 
correction to the gluon propagator with the blob denoting gauge fields/gauginos/scalars/ghosts 
propagating along the loop. Their color factors equal C F N and they directly contribute to 
vj( 2 \ To preserve supersymmetry, we evaluate these diagrams within the dimensional reduction 
scheme. The two-loop correction is given by the sum over the individual diagrams shown 
in Fig. [2] plus crossing symmetric diagrams obtained either by cyclic permutations of the cusp 
points (1,2, ... ,n) i— > (n, 1, . . . , n — 1), or by flips (1, 2, . . . , n) i— > (n, n — 1, . . . , 1) with n — 6. 

To compute the diagrams shown in Fig. [21 we employ the technique developed in Refs. [381 
[19]. The calculation goes along the same lines as for the rectangular and pentagonal Wilson 
loops [TTI [20l |2T] and the result can be summarized as follows. It is convenient to expand the 
contribution of each diagram in powers of 1/e and separate the UV divergent and finite parts as 

- (2) = E \\ ih A - 1 + ^ + 7* A - 2 + I*-}) E ^vf + 4 Q) } + 0(e) , (37) 

5 Up to three loops, the color factors in front of on the right-hand side of (|36|) have the form C F N n and the 
only source of non-planar corrections is the subleading term in the expression for the Casimir C F = N/2 — l/(2N). 
Starting from four loops, the color factors are not expressible in terms of simple Casimir operators and receive 
genuine nonplanar corrections [53] . 
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i — — x 2 

(a) (b) (c) (d) 




(q) (r) (s) (t) 




Figure 2: The maximally non-Abelian Feynman diagrams of different topology contributing to F§ 
The double lines depict the integration contour Cq, the dashed lines the gluon propagator and the blob 
the onedoop polarization operator. 
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where e = e uv and the scale fj,^ Y is defined in ( 1231) . Here the index a runs over the two- 
loop Feynman diagrams shown in Fig. [2(a)-(u) and the factor 1/2 has been inserted for later 
convenience. For each diagram, the coefficient functions depend on dimensionless ratios 
of distances xf i+2 and xf i+3 (with i = 1, ... ,6 and the periodicity condition i + 6 = i). We 
would like to stress that the contribution of each individual diagram to w^ 2 \ or equivalently the 
functions A^l, are gauge dependent and it is only their sum on the right-hand side of (l3~Tj) that 



remains gauge invariant. We note that many of the coefficient functions can be cross-checked 
against similar contributions to W{C§) and WiC^) [2Q112I], which can be obtained by shrinking 
one and two segments, respectively, of the contour Cq to a point. 

3.1 Divergent part 

Summarizing our results for the UV divergent part of W(C e ), we find that the coefficients A^l, 
A^l and A^l 2 of the poles on the right-hand side of (13 7p are given by: 

• 0(e~ 4 ) terms only come from the two Feynman diagrams shown in Figs. 12(a) and[2(r). 

A ( i = -— , A {T \ = — (38) 

-4 lg . -4 lg V J 

• 0(e~ 3 ) terms only come from the two Feynman diagrams shown in Figs. E(r) and[2(s). 

= -i , A% = 1 (39) 

• 0(e~ 2 ) terms only come from the Feynman diagrams shown in Figs. 12(a), (e), (o), (r) and 
(s). 



All these contributions match the corresponding ones for W^Ci) and W (Cj 



5 J 



• 0{e x ) terms come from the Feynman diagrams shown in Figs. E(a)-T5(e) J2(m) ~E(p) and 
0i) 0u). 

The expressions for A_ l7 A_\, A^°\ and A_ t are the same as for W(C$) and W{C$), while the 
remaining coefficients are given by complicated functions of the distances x\ i+2 and x? i+3 . 
To determine them we apply the 'subtraction procedure' described in [21]. As we will see shortly, 
these functions cancel against each other in the sum of all diagrams leading to 

£43 = ^3. (4i) 

a 

Substituting the relations (138]) - (14T1) into (1371) we find that the divergent part of is given by 

{ ^ + e-^Ca} E("^ 2 + ^(e ) . (42) 
^ ' t=i 
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As we have shown in [21], a similar relation holds for an arbitrary n— gon light-like Wilson loop 
W(C n ) (with n > 4). 

In the planar limit, we take into account that Cp = N/2 and rewrite (1361) as 

In W(C 6 ) = aw w + 2aV 2) + 0(a 3 ) . (43) 

Combining it with fj42|) we see that the divergent part of IiiW(Cq) is of the expected form (T2T1) . 

Let us now return to (j4"Tl) and specify the contribution of individual diagrams. We introduce, 
following [21] , the auxiliary Feynman diagrams shown in Fig. [3j These diagrams involve a new 
fake 'interaction vertex' for three gluons. Notice that one of the gluons in Fig. [3] is attached 
to a corner of the hexagon, while the positions of the two remaining gluons are integrated over 
one adjacent segment and one non- adjacent segment. By definition [21], the Feynman integrals 
associated with the two diagrams shown in Fig. [3] are 

4uL = t:9 4 CfN(xi2 ■ x 56 ) I drx f dr 5 J(xi, x 2 + nx 12 , x 6 + r 5 x 56 ) , (44) 
2 Jo Jo 

^ = 7;9 4 CfN(x 61 ■ x 3A ) / dr 3 I dr§J(xi + nx 61 , x x ,x 4 + t 3 x u ) , 

* Jo Jo 

4ux = t:9 4 C f N(x 12 ■ x 34 ) / dr x i dr 3 J {x x ,x 2 + nx 12 , x 4 + 

* Jo Jo 



where we have used the parameterization (Tl8|) of the integration contour Cq and have introduced 
the notation for the auxiliary function [9j [21] 



J( Zl , z 2 , z 3 ) = -t(n 2 uv )- e J d A ~ 2e z G(z - Zl )G(z - z 2 )G(z - z 3 ) , (45) 



with the gluon propagator G(x) defined in flA.3j) . J(z\, z 2 , z 3 ) is a symmetric function of the 



three points zf (with i = 1,2, 3) in Minkowski space-time. 

We observe that on the right-hand side of ( 1441) two of the points, say z 2 and z 3 , are separated 
by a light-like interval, z| 3 = (z 2 — z 3 ) 2 = 0. In this case, J(zi, z 2 , z 3 ) develops a simple pole 
in e (see Eq. f ]A.8j) for the explicit expression). It originates from the integration in fj45|) over 



z^ approaching the light-like direction defined by the vector (z 2 — 23)^, so that the distances 
{z — z 2 ) 2 and (z — z 3 ) 2 vanish simultaneously and the two propagators on the right-hand side 
of fj4"5l) become singular. As a consequence, the integrals liux have a simple pole in e. Then, 
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we add together the crossing symmetric diagrams of the same topology as in Fig. [3J and expand 
their contributions in e similar to (137)) . 



/ 2 \ 2 [" i 6 

J<£ + (cross-symmetry) = C F iV -M^ £ 



2 \ 
2 /\jvJ 



2e + M^+0(e) 



(46) 



with a = {a,b,c}. Here M_" and are complicated functions of the distances x^ fc . We will 
not need their explicit form for our purposes. 

The crucial observation made in [21] is that the coefficient functions A^*[ corresponding to 
the various diagrams shown in Fig. [2] can be expressed in terms of the three functions M^f , 
M^l and Mi(. In Appendix A we illustrate the underlying mechanism by presenting a detailed 
calculation of the diagram shown in Fig. EJ^c). Repeating the similar analysis for the remaining 
diagrams, we obtained the following expressions for the coefficient functions A_\: 



A 



A 



A 



(a) 
1 

(d) 
-1 

to 



24 



A^l 



2M 



(a) 



•1 ■ 



.4 



(c) 



2M 



(b) 



-1 5 



(47) 



2M 



(c) 

-1 5 



.4 



(e) 



A { ™ ] + A 



(t) _ _ M to _ M (c) ^ 



.4 



to 



-M 



to 



A 



(o) 



.4 



(p) 
-i 



. M to _ , 



^w + A {s i = — c 3 . 

-l -l 24 s^ 



The functions M^, and can be expressed in terms of polylogarithms of degree 3. 

Assigning degree 3 to the constant £3, we can say that the expressions appearing in ( 1471) have 
the same degree of 'transcendentality' [§. Adding them together we find that the functions , 
M^l and M_\ cancel in the sum of all diagrams and we recover the relation (EH). 

To summarize, we have demonstrated by an explicit two-loop calculation that the divergent 
part of the hexagon Wilson loop has the expected form (121)) with the cusp and collinear anomalous 
dimensions given by (1T2)) and (122)) . respectively. 



rto 



r(c) 



3.2 Finite part 

The finite part of W(Cq) receives contributions from all Feynman diagrams shown in Fig. [2j 
Substituting (1371) into (136)) and comparing with (1201) . we find the two-loop contribution to Fg WL ^ = 

n P p( WL ) 

^ L) = 2 £ 4°. (48) 

a=a,...,u 

where the overall factor 2 comes from (g 2 / (Ait 2 )) 2 C fN = 2a 2 at large A/". Based on our analysis 
of the rectangular and pentagon Wilson loops [20], we expect that the sum on the right-hand 
side of this relation should have degree of transcendentality 4. 

As was already explained, the contribution of the diagrams shown in Figs. 12(a), (e), (r) and 
(s) to the finite part of In W(Cq) can be obtained from the similar expressions for the rectangular 



3 Here we follow the nomenclature which has become common in the literature after the proposal of 
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Wilson loop [2D], 



4 2880 ' v ; 



4 2880 ' 



.(e) _ 71-2 ^,^2 fgj+M-A 6 19 4 



4 e) = ^Vln 2 f 
96 ^ I 



<i-2 I 4 720 



-7T 



where the combinatorial factor | accounts for the difference in the number of diagrams of the 
same topology between the hexagonal and rectangular Wilson loops. Moreover, the contribution 
of the diagrams shown in Fig. [2](h), (k) and (1) factorizes into a product of two one- loop integrals 
and can easily be calculated. 

Let us separate the remaining diagrams into two groups according to their behavior as e — > 0: 

• The diagrams shown in Figs. [2(f), (g), (i), (j), (q) are finite; 

• The diagrams shown in Figs. E^b), (c), (d), (m), (n), (p), (t), (u) have a simple pole and 
the diagram in Fig. E(o) has a double pole. 

For the second group of diagrams we identify their finite part by employing the 'subtraction 
procedure' [21] ■ This amounts to subtracting from these diagrams the auxiliary diagrams shown 
in Fig. [3] with the appropriate weights defined by the coefficients in front of the functions M_^^ 
in (j47p . In this way, we compensate the simple poles in the above mentioned diagrams and 
render their contribution finite (with the only exception of the diagram in Fig. [2](o), where an 
additional subtraction of the double pole defined in pOl is required). Since the auxiliary diagrams 
also generate a finite part M^' h,c \ Eq. (T4*6|) . the subtractions will modify the finite parts of the 
individual diagrams, — > Aq*\ 

= 4 b) - 2M< a) , A^ = A^ ) -2Ml ) b \ A.f = Af ] - 2M^\ (50) 

^(n)+(u) = A (n) + A {u) + M (b) ? £(o) = A (o) + M (a) ? ftp) = ftp) + M 0>) + ^(c) ? 

j£m)+(t) = ^(m) + A (t) + M (a) + ^c) 

This, however, does not affect the total sum of diagrams, 

£4 a) = £^ a) - ( 51 ) 

a a 

By construction, the subtracted diagrams are free from UV divergences and, therefore, can be 
directly evaluated in D = 4 dimensions. In this way, we find that, 

A^ )+(u) =lj m)+(t) =0. (52) 
In addition, we obtain a representation for the remaining functions 

^(b,c,d,o,p) A (f,g,i,j,q) 

111 

the form of convergent multiple integrals. Being combined with the explicit expressions for the 
remaining functions 7 4^ a ' h ' e,k ' 1 ' I> ) ; they determine the two- loop correction to the finite part of the 
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hexagon Wilson loop ( )48l) . Their explicit expressions are lengthy and we do not present them 
here to save space. As an example, we give the detailed calculation of in the Appendix. 

It is straightforward to evaluate numerically the sum on the right-hand side of (1481) for a given 
set of kinematical invariants, 

7^r2 2 2 2 2 2 2 2 2l /rn\ 
i\ — \x 13 , x 14 , x 15 , x 24 , x 25 , x 35 , x 26 , x 36 , x 46 j , l^Jj 

defined as the distances between the vertices of the hexagon C e . We should take into account, 
however, that in D = 4 dimensions there exists a relationship between the distances (1531) . such 
that only eight of them are independent. This relation reflects the fact that the five four- 
dimensional vectors pf — xf — (with % = 1, . . . , 5) are linearly dependent, and hence their 
Gram determinant vanishes 



G[K) = det \\(p t - Pj )\\ = , (/../ 1 •.). (54) 

Using pi = Xi — the entries of the matrix can be written as linear combinations of the 
distances flMD- 

We would like to note that the sum on the right-hand side of (J48]) defines a function of xjj 
even for configurations K which do not satisfy (|54|) . This function can be viewed as a particular 
continuation of the finite part of the hexagon Wilson loop off the hypersurface defined by (|54|) . 



4 Duality relation for the six-gluon MHV amplitude 

In the parallel publication [13j, an impressive two-loop calculation of the six-gluon MHV planar 
amplitude in M = 4 SYM theory has been performed. To verify the duality relation (I35p we shall 
compare our results for the finite part of the hexagon Wilson loop with the numerical results 
from Ref. [13]. 



4.1 The hexagon Wilson loop versus the BDS ansatz 

Let us first consider the relation between the hexagon Wilson loop and the BDS ansatz, Eq. fl3~4"|) . 
To two-loop accuracy it takes the form 

R w ( Ul , u 2 , t^) = F$ L) - Fg DS) , (55) 

where i?w an d -^6 2 DS ^ denote the two-loop contributions to the remainder function §5§ and to the 

BDS ansatz (fl3j) . respectively. We recall that the functions and satisfy the Ward 

identity ff24l and i?w is a function of the three conformal cross-ratios (|25|) only. 

The simplest way to check relation fl55l) is to evaluate the difference F^ L ^ — F^ DS ^ for two 
different kinematical configurations K and K', related to each other by a conformal SO(2,4) 
transformations of the coordinates (with i = 1, . . . , 6). Since F^ h ^ and F^^^ are dimension- 
less functions of the distances they are automatically invariant under translations, Lorentz 
rotations and dilatations of the coordinates . The only non-trivial transformations are the 
special conformal transformations (boosts), which are combinations of an inversion, a translation 
and another inversion. 
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Let us start with the kinematical configuration K = K(xfj), Eq. ( l53i) . and perform an inver- 
sion of the coordinates, xf — > jx\ , to define the new configuration 



K' = K (4/(x^ 2 )) • (56) 

Since the variables u ( 1251) are invariant under such transformations, w a [-^] = u a [K'], the difference 
~ -^6-2 DS ' > should also be invariant. 
As an example, let us consider six light-like four- dimensional vectors p^, 



p 1 = (l,l,0,0), p 2 = (-l,p,p,0), p 3 = (l,-p,p,0), (57) 

P4= (-1,-1,0,0), p 5 = (1, ~P, ~P, 0) , Pe = (-1,P, ~P, 0) , 

with p = l/v2 and J^pf = 0. These vectors define the external momenta of the gluons in the 
six-gluon amplitude. Applying the duality relation pj = — Xi + x, we evaluate the corresponding 
distances ( 1531) . 

: x 14 = x 15 = x 24 = x 25 = —2 , rc 36 = —2 — 2\/2 , 

X 13 = X 35 = X 26 = X 4S = ~2 ~ V^2, (58) 

and the conformal cross-ratios (1251) . 

Mi = u 3 = \ + fv 7 ^, u 2 = l. (59) 



By construction, this kinematical configuration satisfies the Gram determinant constraint (J|2 
To define the conformal transformations fl56l) . we choose an arbitrary reference four- vector x± = 
(X®, x\, x\ , xf) and reconstruct the remaining a;— vectors according to = — pf (due to 
translation invariance, the Wilson loop does not depend on the choice of x±). Then, relation 
f l53|) implies that the function R w evaluated for the kinematical configuration floT)!) should be the 
same as for the original configuration K. 

As an example, we choose x^ = (1, 1, 1, 1) and apply the conformal transformation (156|) to 
defined in (1581) to obtain the new kinematical configuration 

K {h) . 2 2 _ 2 2 _ _1 _ 1 /o r 2 _ _1 _ 3 

-fV . X 14 — X 15 — X 24 — X 25 — 2 4 V ^ , X 3 g — 1 4 V ^ , 

x 13 — 2 V ' x 35 — 2 4 VZ , x 26 — 4 8 V ' 

x 2 6 = -|-ix/2. (60) 

The results of our numerical tests are summarized in Table [TJ They clearly show that F^ h ^ 

and F 6 ( ; f S) vary under conformal transformations whereas their difference R w = F^ L) - F<g DS) 
stays invariant. 

We recall that in four dimensions the kinematical invariants ( 1531) have to verify the Gram 
determinant constraint ( 15^1) . This relation G[K] = is invariant under the conformal trans- 
formations (I56I) . simply because the conformal boosts map six light-like vectors pf into another 
set of light-like vectors. There exist, however, certain kinematical configurations K' for which 
Ui[K] = Ui[K'] but G[i^'] 7^ 0. Since the difference function i? w only depends on the u— variables, 
its value should be insensitive to the Gram determinant condition 0. For example, consider the 



7 We recall that the functions entering (|55p can be defined for configurations K' satisfying G[K' ] ^ 0. 
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following kinematical configuration 



** : x 14 — x 15 — 1 , x 2 4 — a^25 — 2 , x 36 — 2 2\/2 , 

x? 3 = -1-1/^2, 4; = 4, =4, = -2-^2. (61) 

The corresponding conformal cross-ratios (12"5|) are given by J52J, but G[AT< C )] ^ 0. We verified 
numerically that i? w [^ (a) ] = #w[^ (b) ] = #w[^ (c) ] with accuracy < 1(T 5 (see Table!]). This 
observation allows us to study the function Rw(ui, u 2 , 113) without any reference to the Gram 
determinant condition. 



Kinematical point 


F (WL) 
r 6;2 


p(BDS) 
r 6;2 






-5.014825 


-14.294864 


9.280039 


k o>) 


-6.414907 


-15.694947 


9.280040 


K (o) 


-5.714868 


-14.994906 


9.280038 



Table 1: Two-loop contributions to the hexagon Wilson loop, F 6 . 2 , to the BDS ansatz, F 6 . 2 , 
and to their difference, R^, evaluated for three kinematical configurations (158]) . (160]) and flST]) 
corresponding to the same values of u±, u 2 and u 3 , Eq. (159]) . 



The properties of the remainder function (155]) were studied in our previous publication 
Ref. |24j. We found that R-w(ui,u 2 ,u 3 ) is a completely symmetric function of Mj. For u\ = 0, 
u 2 = u and u 3 = 1 — u its value is independent of u, 

i?w(0, u, 1 — u) = cw • (62) 

As explained in |24j , this property ensures that the hexagon Wilson loop has the same behavior 
as the gluon amplitude in the collinear limit. For the Wilson loop this corresponds to flattening 
one of the cusp points. This lends additional support to the duality relation between the Wilson 
loops and the scattering amplitudes. We found numerically that the value of cw is given by 

c w = 12.1756. (63) 

This number has a much smaller absolute accuracy (~ 10 -3 ) compared to the numbers in Table 1 
because its evaluation requires taking differences of quantities that diverge in the collinear limit. 



4.2 The hexagon Wilson loop versus the six- gluon MHV amplitude 

We now compare numerically the function i?w defined in (155]) with the analogous quantity defined 
for the MHV amplitude, 

R A = 4r V) - F§° S) . (64) 

In order to test the duality relation (135]) . we have to show that for general kinematical configu- 
rations K defined in equation (153]) . 

R W [K] - R A [K] = c w . (65) 
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Since i?w is a function of ui, U2, M3 only, this relation would imply that so is Ra- 

To get rid of the constant &w, which, as we have seen, has lower numerical precision than the 
evaluations for generic kinematics, we subtract from (1651) the same relation evaluated for some 
reference kinematical configuration K^°\ 

R A [K] - R A [K<®] = R W [K] - R W [K®] . (66) 

The numerical tests of this relation for different kinematical configurations are summarized in 
Table d 



Kinematical point 


(ui,u 2 ,u 3 ) 


p d(0) 


Ka — ri A 




(1/4, 1/4, 1/4) 


< lO" 5 


-0.018 ±0.023 


K <2) 


(0.547253, 0.203822, 0.88127) 


-2.75533 


-2.753 ±0.015 


K (B) 


(28/17, 16/5, 112/85) 


-4.74460 


-4.7445 ± 0.0075 


AT(4) 


(1/9,1/9,1/9) 


4.09138 


4.12 ±0.10 




(4/81,4/81,4/81) 


9.72553 


10.00 ±0.50 



Table 2: Comparison of the deviation from the BDS ansatz of the Wilson loop, Rw, an d of 
the six-gluon amplitude, Ra, evaluated for the kinematical configurations (ISTI) . Here, R$ = 
i?w(l/4, 1/4, 1/4) = 13.26530 and R ( A ] = i2 A (l/4, 1/4, 1/4) = 1.0937 ± 0.0057 denote the same 
quantities evaluated at the reference kinematical point K^°\ The numerical results for Ra and 
R A are taken from Ref. [13]. 



: 


x 2 

X i,i+2 


~~ 1 J X i,i+3 ~ 


-2; 










x 13 = 


-0.7236200 , 


x 24 


= -0.9213500 


) ^35 


= -0.2723200, 






x 15 = 


-0.4235500 , 


x 26 


= -0.3218573 


, x 1A 


= -2.1486192, 


X 25 = 




x 2 - 

x 36 ~~ 


-0.4825841 ; 












: 


x 2 - 


-0.3223100, 


x 2 

x 24 


= -0.2323220 


x 2 

' x 35 


= -0.5238300 , 


x 2 - 




x 15 = 


-0.5323200, 


X 26 


= -0.9237600 




= -0.7322000 , 


X 25 = 




X 36 = 


-0.6626116; 












K®: 


T 2 

x i,i+2 


= —1 , a?i4 




1/2, 4, = - 


5/8, 


x 2 36 = -17/14; 




: 


x 2 


~~ 1 5 X i,i+3 ~ 


-3; 








: 


x 2 


= ~ 1 ' ^4,4+3 = 


-9/2. 









-0.3582300, 
-0.7264904 , 

-0.8237640 , 
-0.8286700 , 



(67) 

Among the six configurations in equation (ISTI) only the first four verify the Gram determinant 
condition (|5~4"|) . Also, the configurations and are related to each other by a conformal 
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transformation. This explains why the first entry in the third column of Table [2] is almost zero, 
and also reflects the high precision of the numerical evaluation for the Wilson loop. Comparing 
the numerical values for the six-gluon amplitude and the hexagon Wilson loop, we observe that 
their finite parts coincide within the error bars. Therefore, we conclude that the duality relation 
(j3D is satisfied, at least to two loops. 

5 Conclusions 

The results presented in this paper and in the parallel work [13J show that the conjectured Wilson 
loop/MHV amplitude duality holds at two loops for six cusp points/gluons, extending previous 
checks in a highly non-trivial way. We believe the evidence presented here indicates very strongly 
that the duality holds to all orders in the coupling. As pointed out in [23] , this necessarily implies 
the failure of the BDS ansatz for six gluons and presumably beyond. Furthermore, the fact that 
the Wilson loop and the MHV amplitude both differ from the BDS ansatz by the same non-trivial 
function R w (ui, U2, "U3) suggests that there is much more to the proposed duality than a dual 
conformal symmetry of the gluon amplitude. An issue which needs to be addressed in the near 
future is the analytical evaluation of the remainder function. It could give us clues about hidden 
symmetries of the problem and hints for generalizing it to higher loops, as well as to strong 
coupling. 

While at strong coupling the duality naturally emerges from the gauge/string duality, at weak 
coupling we are still lacking the understanding of its origin. Since we expect that the duality 
applies only to planar amplitudes, this suggests that the best way to address it will not be the 
Lagrangian formulation of Af = 4 SYM. Rather, it is more likely to follow from some effective 
description of planar MHV amplitudes. A possible candidate might be the twistor approach 
proposed in [55] . 

Finally, let us return to the dual conformal symmetry of the gluon amplitudes. Assuming 
that the duality with Wilson loops holds, this symmetry follows from the ordinary conformal 
symmetry of the light-like Wilson loops in Af = 4 SYM. This does not explain, however, the 
origin of the duality itself. We believe that further understanding of the duality relation will 
come from a deeper investigation of the symmetries of gluon amplitudes. One might speculate 
that dual conformal symmetry is just part of a larger set of symmetries yet to be discovered. One 
possible way to uncover such symmetries might be the loop equations [5HE7] both at weak and 
strong coupling [57', '58] . Were these symmetries powerful enough to determine the finite parts of 
the planar amplitudes to all loops, one would conclude that this sector of Af = 4 SYM is exactly 
solvable. 
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A Appendix: Two loop corrections to the hexagon Wil- 
son loop 

In this Appendix we present some details of the Wilson loop calculation. It serves to illustrate 
both the technique employed and to outline the general structure of two-loop corrections to the 
finite part of the hexagon Wilson loop. We will examine one diagram in detail, while for the rest 
we will just summarize the relevant formulae. 

As our detailed example, we consider the Feynman diagram shown in Fig. [5(c). It involves 
two gluon propagators and originates from the expansion of the path-ordered exponential (jTTj) 
to fourth order in gauge fields. The end-points of the propagators are attached to three light-like 
segments [x6,xi], [^1,^2] and [2:3,2:4]. It is convenient to parameterize the points belonging to 
the light-like segment x^{t 3 ) G [2^,2^+1] in the following way 

x^rj) = x» +1 + r 3 a^ +1 , (0 < r 3 < 1) (A.l) 

with Xjj+i = Xj — Xj +l and x 2 j +1 = 0. Then the contribution of the diagram shown in Fig. [2(c) 
to the Wilson loop reads 

I c = (ig) A N~ l tx[t a t h t a t h ] f dnx^ f dr 6 2$ G mw (x 12 (l - n) + x 61 r 6 ) (A.2) 

Jo Jo 

dr 3 x%\ I dr' G Xq{ G^ U6 (x 13 + 234(1 - r 3 ) + x 6 iTq) 





where the condition r' % < r 6 is due to the ordering of the gluons along the path and G^ u (x — y) 
is the gluon propagator. In our calculation we employ the Feynman gauge and use dimensional 
regularization with D = 4 — 2e and e > 0, so that the gluon propagator is given by 

G^ix) =g^G(x), G(x) = ~ e) (-x 2 + »0)- 1+e (/x 2 7r) e . (A.3) 

We would like to stress that the contribution of each individual Feynman diagram to the Wilson 
loop is gauge dependent whereas the sum of all two-loop diagrams shown in Fig. [2] is gauge 
invariant. As a first step, we evaluate the color factor of I c 

N- 1 ti[t a t b t a t b ] = C F (C F - §JV) , (A.4) 

with Cp = t a t a = (N 2 — 1)/(2N) being the quadratic Casimir of SU(N) in the fundamental 
representation, and verify that it involves the maximally non-Abelian part ~ CpN. As explained 
in Sect. 3, it is only this part of I c which contributes to hxW(Ce). Then, we make use of the 
identities 2^ - +1 = to simplify [212(1 — Ti) + xqiTq] 2 = 2(si 2 26i)(l — t\)tq and similarly for the 
second propagator in (IA.2I) . The integration over n, r 3 and tq can easily be performed leading 
to 



%a = ( ^3 ) CpN^^i^Ti-^Y (A.5) 



X 



/ dt (t e - 1)— 2 2 \, T 7Z3, Z2~\r [{~ x iJ - xl 6 t) - (—x 13 t - x 2 m t) ] 
Jo yx m x 3e )L -1- (x 14 x 13 )L 
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with t = Tg and t = 1 — t. Here the subscript indicates that we retained only the maximal 
non-Abelian part of I c . We observe that I^ A has a single pole as e — > 0. It comes from the 
integration region n — > 1. 

To separate the divergent and finite parts of J^a we use the subtraction procedure described 
in detail in Ref. [21]. We introduce the auxiliary Feynman diagram shown in Fig. [3]^b), 

J iui = t;9 4 C f N(x u ■ x 61 ) / dr 3 I dTQ J(x x + t 6 x 61 ,x 1 ,x 4 + t 3 x u ) , (A. 6) 

* Jo Jo 

where J(zi, z 2 , z 3 ) stands for three propagators joined at the same point z which is integrated 
out , 

J{z 1 ,z 2 ,z 3 ) = -i(fi 2 )~ £ I d D zG{z-z l )G{z-z 2 )G{z-z 3 ). (A.7) 



The relation (1A.6I) involves this function evaluated for z\ = Xi + r 6 x 6 i, z 2 = X\ and z 3 = x 4 + t 3 x 34) 
so that the points Z\ and z 2 are separated by a light-like interval. For (z\ — z 2 ) 2 = the calculation 
of J(zi, z 2 , z 3 ) yields [9] 



2 . 2£ r(l-2 £ ) f 1 dr(rf) 



J( Zl ,z 2 ,z 3 ) = (7r/i 2 ) 2£ \ / r7 , 211 _ 2£ (A.8) 

with f = 1 — r. We substitute this relation into flA.6j) . change the integration variable via t = r e r 
and expand liux in powers of e to find after some algebra, 

The integral on the right-hand side of this relation can be expressed in terms of polylogarithms 
of degree 3. Symmetrizing Zaux with respect to cyclic permutations of indices, £—►« + !, and 



flips, i — ► 7 — i, we arrive at (j46j). 

Let us now compare the expressions for the two diagrams, Eqs. AA.5D and flA.9[) . We observe 
that both iIsmx and I§ A contain a single pole 1/e, but it disappears in the combination I§j A — 2I^ x . 
This suggests writing 7§a as 

/£ A = 2/i£ + / fi c n + 0(e) (A.10) 
with Jg n finite as e — > 0. Combining (1A.5[) and (1A.9j) we find 

4 C n = -5 f fV) ^ - ^ - + 4) (A.H) 
8V47TV V^26 X 26 X 26 X 26 J 







1 dt ln(Xi(t)/X 2 (t)) 

Xt(t)-X 2 (t) 



llntln(X 1 (t)X 2 (t)/t 2 )+ / ^ln(l-y) 

2 Jt y 



where 

X,(t) = ^(1 - t) + ^ , X 2 (t) = ^(1 - t) + ^ . (A.12) 

x 26 x 26 x 26 X 2G 

Taken together, the relations (1A.10I) . ( 1A.11I) and ( 1A.9I) define the contribution of the diagram 
shown in Fig. EIc) to \to.W(Cq). Symmetrization of Jg n with respect to cyclic permutations of 
indices and flips yields the expression for A ( c) defined in fl50]). 
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We observe that both the divergent and finite parts of I^ A , Eq. ( 1A.10I) . are complicated 
dimensionless functions of four ratios x^q/x^q, x\§jx\^ x\Jx\§ and x\ z jx\§. As such, they are not 
invariant under conformal transformations of the coordinates xf (with i = 1, . . . , 6). Moreover, 
close examination of ij^L and shows that the integrals entering the relations flA.91) and flA.111) 
have transcendentality 3 and 4, respectively. Assigning one unit of the transcendentality to 1/e 
we find that has degree of transcendentality 4. 

We will now summarize the remaining non-trivial contributions to the finite part of lnW^g. 
Each one must be included in all of its inequivalent orientations. Diagrams (a), (r) and (s) are 
identical to those appearing in the calculation of the rectangular Wilson loop and they contribute 
only a constant contribution to the finite part. Their explicit expressions can be found in p^9| [20] . 
The combinations (m) + (t) and (n) + (u) also do not contribute to the finite part, as described in 
[2"T] . For the remaining diagrams, some of the contributions (diagrams shown in Fig. 2 (h),(k),(l)) 
factorize into products of finite one-loop integrals given by 



J 



dtjdtk (xjj + \Xk,k+i) 

I) ( x jk ~ x j,j+ltj + x k.k+ltk) 2 



1 



dt- 



x 



•2 ™ 2 \ / ,-v»2 ™ 2 

j'+l,fc+l ^j^k+l) K^j+l^k ^jh 







X 



J+l,fe+l 



X :.. , i)t + {x 2 - +l k 



c 2 



= ln 



X 



j+l,k+l 



t 



X j+l,kt 



X 



j,k+l 



t+x ikt 



(A.13) 



For the remaining expressions we will often make use of the following shorthand notation for 
certain ratios of the variables x?-, 



a: 



2 2 2 2 2 

_26 Q _ ^36 _ -HA s. _ Z46 _ x 15 

__2 ' P __9 J ' _„2 ' __2 ' '/ __2 " 



X 



13 



X 



13 



X 



13 



X 



13 



X 



13 



(A.14) 



Sometimes it is also convenient to use the notation (recall that Pi 



Xi+l) 



Sij 



^Pi ' Pj 2Xj j j-)-i • Xjfj'+l 



2,2 2 

rf _|_ r*> rr> 



X 2 

x i+l,j+l- 



(A.15) 



Further, we will often use the notation y = 1 — y for some variable y. We will never use the bar 
to refer to complex conjugation. Using this notation, the remaining non-trivial contributions are: 

Diagram (b) 



U 9 



'fin 



2 \ 2 



C F N 



dt 



8 V47T 2 / Jo t + S 12 /S 2 6 

where X(t) = I /(at) + 1 + s 2 §/si & . 



iln 2 X(t)lnt + lnX(t) — \ny 
2 Jt V 



(A.16) 



Diagram (d) 



J fin 



l 



hi 2/ 



(A.17) 



where X\ 



x 
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t)/(a& - a4 - 2&) and X 2 = (x 2 41 t + x 2 m t)/(x 



x 



1(5 



x 5l)- 
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Diagram (e) 



fin 1 4tt 2 



g2 ^C F {- l -N) 



C2 j n 2 ( x 15 



19 ^ 



X2 6 J 1440 



(A.18) 



Diagram (f) 



v] [5 - A - a - 1 - 1] 



(A.19) 



x In 



A(A + 1 - f) + a (A -t))i r - t) - t(a + A) 



a(X — u) 

where A = (3 — a — 1 and /x = 1 + 5 — 7 — (3. 



In 



/x(A + a + 1) - t(a + A) 



Diagram (g) 



P 



■ J ' ; C F (-\N) f dt [" t/r 



in 



2\2 



1 



(l Jo lV -\-i][5 + t] 
A(A+ 1 - v) + a(A - v)- 



(A.20) 



x In 



a(A — f) 



In 



p5 + (p + 77)*. 



In the expression for P the following notation has been used: A = (3 — a — 1 and p = 7 — 77 — 5. 
Diagram (h) 



(4vr 



2^2 



C F (-\N)J M J, 
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(A.21) 



where is given by flA.13j) . 
Diagram (i) 



57T 



2^2 



C F (-\N) [ dv [ dt—^ — 
2 Jo Jo v + Xt- 



(A.22) 



x In 



t(v + A) +a(\ - v) 
a(A — v) 



In 



t> (t — a) + tX + Aa 



In the expression for P the following notation has been used: A = 1 + a — (3 and A = 1 — A. 
Diagram (j) 



P 



>7T 



2\2 



C F {-\N) dv dt 



jo 



[6-p-v][l-P-t] 



(A.23) 



x In 



t(5-{3-v)+Pp + v(l-(3) 
Pn + v(l-p) 



In 



/i(l - t) - t(/3 - (5) 
Vv{\ - 13 - t) + pn - t{P - 5)\ 



In the expression for P the following notation has been used: /x = 1 — /3 — 7 + 5. 
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Diagram (k) 



Diagram (1) 



(4tt 2 ) 



; C F {-\N)J U J : 



36 ■ 



I 



;C F (-\N)J 15 J 26 . 



( 47r 2)2— V 2 

Diagram (o) 

There are several contributions to the finite part coming from this diagram. We have 

jo _ TO | to | ro 
-'fin — -'fin,! + i fin,2 + i fin,3- 

The first is 



[(ati + at 2 )v + bt 1 t 2 v\ 1 + (a 



The second part is 



4tt 2 7 16 



#i = -2(1 -a -b) [ dt^dxdv--^^ 

Jo vxtib- 

,i x 

K 2 = (1 — a — b) I dt\dt 3 dxdv 

Jo xB - 



v at 3 + bt\ 
xt 3 a 
B 



x 



x 



vbti + vxC vat 3 + vxC 
A 



K 3 = (1 — a — b) / dt\dt 3 dxdv< — 
Jo I x 



Bv + vxC Av + vxC 



xA 1 1 

+ 



xB -xABv + vxC Bv + vC 



1 

+- 

x 



xB 1 1 

+ 



xB - xAAv + vxC Av + vC 



Here we use the notation, 

a = 1/(3 , b = a/ (3 , A = at 3 , B — bt 1 , C = at 3 ii + &M 3 + M 3 . 
Finally the third part is 

•§26 



fin ' 3 32 V4tt 2 J 



2 \ 2 «1 

CpN I dt 3 



S\q + S26^3 



In 



Sl6 + ( s 12 + 526)^3 , 2. , / <_ \ 

In t 3 + (s 12 ±+ s 16 ). 

Sl2^3 
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Diagram (p) 

There are several contributions to the finite part coming from this diagram. We have 

= 'L,l + 11,2 + 'L,3 • (A.32) 

The first part is 

^ = -K£) W t ix [ <, + + *4 ) ■ (A 33) 

with 

Z12 =P2+ P3V2 , Z23 = X31 ~ P3V2 ~ P&v 3 , z 31 =pi+ p 6 v 3 ■ 
Similarly the second is 

= 4 (io w **** i 1 ix r fo+1)( ^ + ^^4) ■ < AM) 

with 

Zl2 = P2V1 + p 3 , Z 2 3 = X41 - P6V3 , Z31 =Pl+ P2V1 + P(jV 3 . 

The final contribution is 

ifins = 7:l T-z) C F N [ dv 1 dv 2 dv 3 [ dx [ - — % ^ 31 + f ^ 2S — 5—, (A. 35) 

fin < 3 8V4ttV J Jo Jo l+yiyxxz^ + xz^ + xzl] 2 ' 

where 

Zl2 = P2Vl + P3V2 , Z 23 = X 31 - p 3 V 2 - p 6 V 3 , Z 31 = pi + J5 2 Wi + p 6 W 3 = -X31 - p 2 ^l + P6^3 

and 

V 23 = -2(p 2 ■ P3){p% ■ X31) + 4(p 3 • Pe)(P2 ■ X31) ~ 2(jo 2 • P3)(P3 ■ Pe)v2 - 4(j» 3 • Pg)(P2 ■ Pe)v3 
V31 = -2(p 2 • p 3 )(pi ■ p 6 ) + 4(p! • p 3 )(p 2 ■ p 6 ) + 2(p 2 • p 3 )(p 2 ■ p e )vi + 4(p 2 • p 6 )(p 3 ■ p 6 )v 3 

Diagram (q) 

g A C F N f 1 



I q = 



1-2 



1287T 4 / dtldt ^ dt ^( tl + *3 + *5 - 1) y n<iu*M3^l3 + *3^35 + *5M5l]' 

x T 3 («9 Z1 , 9 23 , <9 Z J + t 3 t 5 zl 5 + ht x z\ x \ . (A.36) 

Here T 3 is the differential operator, 

r 3 (d zl ,d Z3 , d Z5 ) = -(pi -p 3 )(P5 • (<9 23 - 9 Z1 )) - (p 3 -p 5 )(Pi ■ (d Z5 - d Z3 )) - (p 5 -pi)(p 3 • - <9 Z5 )), 
and the zj k variables are defined by 

z 13 = (Pl^l + P2 + J5 3 V 3 ) 2 , Z 35 = (j9 3 W 3 + p 4 + P5^5) 2 , z\ x = (p 5 V 5 +p 6 + PlVlf . 
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